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SUMMARY 

A  general  steady  throu^-flow  theory  of  nonviscous  fluid  in  turbo- 
machines' of  arbitrary  hub-  and  casing-wall  shapes  with  subsonic  or 
supersonic  velocity  is  presented.  The  theory  is  applicable  to  both 
direct  and  inverse  problems  and  is  derived  primarily  for  use  in  turbo¬ 
machines  having  thin  blades  of  hi^  solidity  with  a  simple  approximate 
correction  factor  for  blade-thickness  effect.  Through  the  use  of  the 
stream  function,  the  continuity  equation  and  the  equation  of  motion  in 
the  radial  direction  are  combined  to  form  a  principal  equation  for  the 
present  problem.  The  principal  equation  contains  some  terms  that  are 
either  prescribed  or  to  be  determined  by  other  equations  defining  the 
problem.  Two  forms  of  the  principal  equation  are  obtained  for  the  two 
main  groups  of  current  compressor  and  turbine  design  in  which  the  vari¬ 
ation  of  tangential  velocity  and  the  variation  of  the  ratio  of  tangen¬ 
tial  to  axial  velocity  throughout  the  blade  region  are  given.  When  the 
tangential  velocity  is  given,  the  principal  equation  is  elliptic  or 
hyperbolic,  depending  on  whether  the  meridional  velocity  is  subsonic  or 
supersonic,  respectively.  When  a  relation  between  the  tangential  and 
the  axial  velocity  is  given,  however,  the  principal  equation  becomes 
hyperbolic  when  the  relative  velocity  is  supersonic.  A  general  method 
of  solution  for  both  the  elliptic  and  the  hyperbolic  case  is  outlined. 
Specific  applications  of  the  theory  to  several  common  types  of  com¬ 
pressor  and  turbine  employing  free-vortex,  symmetrical-velocity-diagram, 
solid-rotation- type,  nontwist ed-blade,  and  radial-blade-element  designs 
are  discussed. 


INTRODUCTION 

With  the  increasing  use  of  velocity  diagrams  other  than  free- 
vortex  type,  low  inlet  hub-tip-radius  ratios,  and  high  velocity  of  flow, 
the  problem  of  three-dimensional  flow  in  axial-flow  turbomachines 
becomes  more  and  more  important.  This  problem  is  treated  by  Traupel, 
Meyer,  and  Marble  (references  1  to  3)  for  incompressible  fluid. 
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CompressilDle  flow  is  treated  in  reference  4,  in  which  methods  for  ^ 
limiting  solutions  of  zero  and  inf inite-hlade  aspect  ratio  are  ohtaine 
and  a  step-hy-step  method,  as  well  as  a  simpler  method  based  on  an 
approximate  knowledge  of  the  shape  of  streamlines,  is  given  for  fini  e- 
hlade  aspect  ratio.  In  reference  5,  Reisner  gives  a  method  of  blade 
design  for  compressible  flow  with  the  shape  of  neither  hub  nor  casing 
wall  specified  in  advance.  The  problem  of  supersonic  flow  in  impellers 
of  a  given  casing  and  blade  shape  is  currently  being  investigated  by 
Arthur  W.  Goldstein  of  the  MCA  Lewis  laboratory. 


The  analysis  made  at  the  MCA  Lewis  laboratory  and  presented 
herein  proposes  a  unified  theory  that  is  applicable  to  both  direct  and 
inverse  problems  for  both  subsonic  and  supersonic  flows  in  compressors 
and  turbines  of  arbitrary  hub-  and  casing-wall  shapes. 


Equations  of  motion  and  energy  for  unsteady  three-dimensional  flow 
of  a  nonvlscous  compressible  fluid  are  expressed  in  terms  of  some  con¬ 
venient  quantities  for  analyzing  flow  in  turbomachines.  Entropy  changes 
due  to  heat  transfer  in  a  cooled  turbine  and  due  to  shock  wave  m 
supersonic  flow  can  easily  be  accommodated.  The  condition  under  which 
irrotational-flow  analysis  is  correct  is  also  discussed. 

The  general  equations  are  then  simplified  for  steady  throu^-flow 
in  turbomachines  having  thin  blades  of  high  solidity. _  It  is  shown 
that,  in  the  direct  problem,  just  enou^  equations  exist  to  detemlne 
all  the  variables;  whereas  in  the  inverse  problem,  after  the  inclusion 
of  the  Integrabillty  condition  for  the  blade  surface,  either  one  vari¬ 
able  or  a  relation  between  several  variables  can  be  specified  by  the 
designer . 

In  the  solution  of  the  problem,  the  continuity  equation  and  the_ 
equation  of  motion  in  the  radial  direction  are  combined  into  a  principal 
equation  through  the  use  of  the  stream  function.  This  equation 
involves  some  terms  that  are  either  given  or  to  be  determined  by  other 
equations  defining  the  present  problem.  Two  forms  of  the  principal  ^ 
equation  are  obtained  for  two  main  groups  of  designs  in  which  the  vari¬ 
ation  of  tangential  velocity  and  the  variation  of  the  ratio  of  tangen¬ 
tial  and  axial  velocities  are  prescribed  by  the  designer.  The  crite- 
rions  of  whether  the  principal  equation  is  elliptic  or  hyperbolic  are 
obtained  for  both  groups. 

A  general  method  of  solving  the  set  of  equations  for  both  the 
direct  and  the  inverse  problems  is  then  described  for  turbomachlnes  of 
arbitrary  hub-  and  casing-wall  shapes  and  for  either  an  elliptic  or 
hyperbolic  principal  equation. 
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*  The  dimensionless  forms  of  the  principal  and  other  equations  for 

some  typical  designs  are  given.  A  simple  approximate  correction  for 
blade  thickness  is  also  given. 

SYMBOLS 

The  following  symbols  are  used  in  this  report: 

a  velocity  of  sound 

m  1 
n^j 

b,f,J,K, 

L,M,W 

C  constant 

Cp^c^  specific  heat  of  gas  at  constant  pressure  and  volume^ 

respectively 

,  ^  differentiation  with  respect  to  time  following  motion  of 

Dt  fluid  particle 

V  derivative  of  q 

F  hlade  force  per  unit  mass  of  fluid 

G  Green’s  function 

y2 

H  total  enthalpy  per  unit  mass  of  fluids  h  + 

h  enthalpy  per  unit  mass  of  fluid 

I  =  h  +  ^-  ^=  H-  CD(V^r) 

k  thennal  conduct ivity 

I  characteristic  curve 


differentiation  coefficients  in  equation  (70)  used  to  multi¬ 
ply  function  value  at  point  Xj  to  give  m*^^  derivative 
at  Xj^  using  polynomial  of  n'^^  degree 

functions  of  r  and  z 


m 


order  of  derivative 
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n 

P 

Q 

q. 

R 

m  i 
nR 

t,6,z 

r, Cp,z 
s 

T 

t 

U 

u 


V 

W 

X 

P 

s. 

r 

r 


degree  of  polynomial 
pressure 

heat  added  to  fluid  particle  along  path  of  motion  per  unit 
mass  per  unit  time 

dependent  variable 

gas  constant 

remainder  term  of  mth  derivative  at  point  Xi  obtained  hy 
using  polynomial  of  n"^^  degree 

cylindrical  coordinates  relative  to  stator 

cylindrical  coordinates  relative  to  rotor 

entropy  per  unit  mass 

temperature 

time 

velocity  of  blade  at  radius  r 

internal  energy  per  unit  mass  with  0°  absolute  as  base 
temperature 

absolute  fluid  velocity 
relative  fluid  velocity 
independent  variable 
Wu 

arc  tan  ^ 

grid  spacing  in  r-  and  z-directions^  respectively 

ratio  of  specific  beats 
average  value  of  X 
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X 


V 
P 

CO 

Subscripts : 
b 
c 
d 


J 

m 

o 

r^u^z 

T 

t 


independent  variables 

dr 

slope  of  characteristic  curve^  — 

n  Wu 

■tan  P  =  ^ 

value  of  X  between  and  x^ 

mass  density 
stream  function 
angular  velocity  of  blade 

trailing  edge  of  guide  vane 
leading  edge  of  rotor 
trailing  edge  of  rotor 
leading  edge  of  stator 
hub 
inlet 

any  station 
meridional 

refers  to  position  where  blade  element  is  radial  or 
radial;  circumferential;  and  axial  components 
total  state 
at  tip  or  casing 


=  0 


Superscripts : 

a;b;  grid  points 


dimensionless  values 
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GENERAL  BASIC  EQUATIONS 

The  three-dimensional  compressible  flow  of  a  nonviscous  gas 
through  a  turbomachine  is  governed  by  the  following  set  of  basic  laws 
of  aerothermodynamics ;  From  the  principle  of  conservation  of  matter^ 
the  equation  of  continuity  is 


or 


D(loge  p) 


(la) 


For  a  blade  rotating  at  a  constant  angular  speed  a>,  Newton's  second 
law  of  motion  gives 

(2) 

The  first  law  of  thermodynamics  may  be  written  as 


Du 

Dt 


+ 


P 


=  Q 


where  u  is  related  to  the  gas  temperature  T  by 


Du  _  OT 
Dt  Dt 


(3) 


(4) 


and  Q  is  given  by  the  following  equation  if  only  conduction  is  con¬ 
sidered: 

Q=p"^y.(kVT)  (5) 

For  the  range  of  gas  temperature  and  pressure  encountered  in  ordinary 
turhomachines ^  p^  and  T  are  accurately  related  hy  the  following 

equation  of  state 

p  =  RpT  (6) 


W 
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Although  the  flow  of  gas  through  the  turbomachine  is  completely 
defined  by  the  preceding  equations  together  with  the  known  variations 
of  Cy  and  k  with  temperature  and  the  given  boundary  and  initial 
conditions ;  in  reference  4  it  is  found  more  convenient  to  express  the 
state  of  the  gas  in  terms  of  entropy^  total  enthalpy^  or  a  quantity  I 
of  the  gas^  in  addition  to  its  velocity  components.  These  quantities 
are  defined  as  follows: 


T  ds  =  du  +  p  d(p"^j 

(7) 

H  =  h  .  ^ 

(8) 

y2 

I  =  h+  -^ - 2  ~  ^  ~ 

(9) 

where 

h  =  u  +  pp  ^ 

(10) 

When  equations  (4),  (6),  (7),  (lO),  and  the  relation, 

R  =  Cp  -  Cv  =  (r-i)c^ 

which  follows  from  equations  (4),  (6),  and  (lO),  are  used, 
obtained 

there  are 

T  ds  =  dh  - 

P 

(11) 

T  ds  =  df£\  -  ^ 

r-1  \pj  p 

(lla) 

‘^(i)  "  7^  p)  -  <3.(loge  p) 

(12) 

<i(|)  =  <i(loge  T)  -  d(loge  p) 

(12a) 

and  the  equation  of  continuity  can  he  written  as 

r-,  7j  .  1  ^(loge  ®  n 

(13) 
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From  equations  (9)  and  (ll)^ 


-  yp  +  i  VW^  -  oo^r  =  VI  -  TVs 

P  2 

With  this  equation  and  the  relation 

™  ^  +  (W.V)W  -  f  +  I  VW2  -  WX(VX¥) 

the  equation  of  motion  (2)  can  he  written  as 

^  _  WX^XW)  +  2^W  =  -VE  +  TVs  0 

An  alternate  form  of  equation  (l4)  that  invplves  the  ^o^icity  of  the 
absolute  motion  is  obtained  as  follows:  Using  the  cylin  ri 
dinate  sys’beui  'wi'bh.  “the  z-axis  parallel  ho  co,  yiel  s 

V  =  W  +  mxr 

yXV  =VXW  +Vx(a)xr) 


VX  (Sxr)  =  (r-V)co  -  (ai-V)r  +  CD(V-r)  -  r(V-a^)  =  2q 


Hence ^ 


V  XV  =VXW  +  2(D 


(15a) 


This  relation  can  also  be  seen  from  the  following  expressions  of  rel  - 
tive  and  absolute  vorticity  expressed  in  terms  of  cylindrical  coo 
dinates  r,cp,z  and  r,e,z,  which  refer  to  the  rotor  and  stator, 

respect ively: 


(Vxw)^  =  7  ^  ^ 


BW-„  Sw„ 


=  Si"  "  ^ 


(VXW)z  =  - 


1  1 
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The  preceding  equations  lead  to  several  important  general  considerations. 

If  the  blade  rows  are  not  placed  too  close  together^  the  pressure 
of  gas  at  a  fixed  point  relative  to  the  blade  can  be  taken  as  constant 
with  respect  to  time.  Consequently^  according  to  equation  (20)^  the 
entropy  and  the  quantity  I  of  the  gas  stay  constant  along  its  rela¬ 
tive  streamlines  for  adiabatic  flow.  The  constancy  of  I  means  that 
the  rate  of  change  of  total  enthalpy  along  the  streamline  is  equal  to 
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the  angular  speed  of  the  hlade  multiplied  hy  the  rate  of  change  of 
angular  momentum  about  the  machine  axis  of  the  gas  particle  along  its 
streamline,  which  is  the  well-known  Euler  turbine  equation  usually 
derived  under  less  general  conditions.  In  a  cooled  turbine  where  the 
heat  transfer  may  be  large,  the  rate  of  change  of  s  and  I  along 
the  streamline  can  be  obtained  by  using  equation  (20).  Again,  for 
steady  relative  flow,  equation  (l4a)  shows  that  the  vanishing  of  abso¬ 
lute  vorticity  requires  both  gradient  I  and  gradient  s  to  vanish 
or  the  difference  between  VI  and  TVs  to  vanish.  When  both  gra  len 
I  and  gradient  s  are  zero  upstream  of  the  blade  row  and  the  flow 
is  adiabatic,  s  is  uniform  in  passing  through  the  blade  row;  p  is 
then  a  function  only  of  p  (according  to  equation  (12)),  and  conse¬ 
quently,  according  to  Kelvin's  circulation  theorem,  the  absolute  vor¬ 
ticity  will  remain  zero  in  passing  throu^  the  blade  row  and  the  flow 
can  then  be  treated  on  the  basis  of  irrotational  absolute  flow. 


For  flow  through  a  stationary  blade  row,  equation  (l4a)  becomes 


VX(V><V)  =  -VH  +  TVs 


(14b) 


which  agrees  with  similar  relations  previously  obtained  by  Vazsonyi 
(reference  6),  and  Hicks,  Guenther,  and  Wasserman  (reference  7).  It 
is  interesting  to  _see  that,  for  relative  flow  in  a  rotating  blade  row, 
VX(VXV)  becomes  WX(VkV)  and  H  becomes  I. 

When  it  is  assumed  that  the  fluid  enters  the  machine  with ^uniform 
H  and  s  and  zero  vorticity,  the  adiabatic  flow  through  the 
guide  vanes  can  be  treated  on  the  basis  of  Irrotational  flow.  When 
the  guide  vanes  impart  a  radial  variation  of  the  tangential  velocity 
of  the  fluid  downstream  of  the  vanes  similar  to  that  in  a  potential 
vortex,  the  circulation  is  constant  along  the  blade  span  and  the  flui 
maintains  a  uniform  H  and  s  and  a  zero  vorticity  going  into  the 
following  rotor  blade  row.  If  the  rotor  blade  row  is  situated  far 
away  from  the  inlet  guide  vanes,  the  fluid  enters  the  rotor  with  a 
uniform  I  in  the  circumferential  direction  as  well  as  ^  the  radi 
direction,  and  the  flow  through  the  rotor  blades  can  again  be  ideated 
on  the  basis  of  zero  absolute  vorticity.  If  the  rotor  is  close  to  th 
guide  vanes,  the  fluid  entering  the  rotor  blades  is  circi^erential  y 
nonuniform,  which  condition  is  balanced  by  the  unsteady  term  of  the  rel 
ative  velocity,  and  the  flow  through  the  rotor  blades  should  theoreti¬ 
cally  be  treated  on  the  basis  of  unsteady  flow  with  zero  absolute 
vorticity.  When  the  guide  vanes  impart  a  radial  variation  of  tangen¬ 
tial  velocity  of  the  fluid  downstream  of  the  vanes  different  from  th 
in  a  free  vortex,  however,  the  circulation  varies  along  the  span  of 
the  guide  vanes,  thereby  shedding  vortices  from  the  trailing  edge  to 
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the  fluids  downstream;  and  the  fluid  enters  the  following  rotor  blades 
with  a  uniform  s  and  H;  a  nonuniform  I;  and  a  nonzero  value  of 
absolute  vorticity.  Consequently;  the  flow  through  the  rotor  blade 
row  cannot  be  treated  on  the  basis  of  zero  absolute  vorticity. 

From  the  preceding  discussion;  the  choice  of  H  or  I  and  s  as 
the  basic  thermodynamic  variables  of  the  gas  is  apparent.  Compressor 
and  turbine  rotors  are  usually  designed  to  add  or  subtract  the  same 
amount  of  energy  radially  to  or  from  the  gas;  hence;  H  is  usually 
radially  constant  throughout  the  machine.  When  the  circumferential 
velocity  of  gas  upstream  of  the  blade  row  is  zero  or  varies  inversely 
with  radius ;  I  is  then  radially  constant  throughout  the  machine  when 
the  heat  transfer  is  zero  or  is  uniform  radially.  These  facts  will  be 
utilized  in  the  following  developments. 

The  continuity  equation  (13)  can  also  be  put  into  a  form  contain¬ 
ing  H  or  I  by  use  of  a  constant  value  of  y*;  that  iS;_  X  for  the 
range  of  temperatures  involved  in  the  process.  By  using  x?  equation 
(l2a)  can  be  written  as 


d(logg  p) 


=  d 


logp 


Integrating  from  the  inlet  total  state  yields 


(21) 


1  1 


where  As*  =  s*-s'*'^  The  continuity  equation  then  takes  the  following 
form: 
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The  solution  of  the  three-dimensional-flow  problem  thus  consists 
of  the  simultaneous  solution  of  equations  (l)  or  (23),  (14),  and  (20) . 

GEKERAL  THROUGH-FLOW  THEORY 

Because  of  the  enormous  mathematical  difficulty  in  solving  the 
preceding  set  of  general  equations,  the  essential  feature  of  the  three- 
dimensional  flow  in  turhomachines  will  he  investigated  with  the  follow¬ 
ing  two  simplifying  conditions:  (l)  The  blade  rows  will  be  assumed  to 
be  placed  so  far  apart  that  the  relative  flow  throu^  any  blade  row  is 
steady.  Under  this  condition,  the  equations  of  continuity,  motion, 
and  energy  in  the  scalar  forms  are : 

Continuity  equation,  from  equation  (l), 

1  3(p¥ji-)  1  a(pwj  SCpWj.) 


or,  from  equation  (23), 


The  three  equations  of  motion,  from  equations  (14), 


W„rS(V,,r)  SW. 


Wj,  5(V^r) 

~  §r  "  Bcp 


^  V Sr 


fj  -  'ST 

,,  /i  Sw^\ 

-  (^r  ^ 


1  _  Sw^ 

r  Sep  Sz;  , 


+  W,,  ^ 


1  Si  T  Ss 
r  r  5^ 


Sz  Sz 


and  the  energy  equation,  from  equation  (20), 

(  Ss  ^u  Ss  Ss  \  Si  ^u  Si  „ 

Q=T(^Wr^  +  —  ^+W2^J=Wr^  +  —  ^+w 


.  Si 
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(2)  The  hlade  is  so  thin  and  the  number  of  blades  is  so  large  or 
the  pitch  is  so  small  compared  with  the  blade  chord  that  the  terms  con¬ 
taining  the  circumferential  variation  of  the  velocity  components  in 
eq.uations  (24)  to  (28)  are  much  smaller  than  other  terms  in  the  same 
equations  and  can  therefore  be  neglected.  The  circumferential  vari¬ 
ation  of  pressure  or  enthalpy  is^  however,  preserved  in  equation  (26) 
by  the  introduction  of  a  blade  force  F,  which  may  be  considered  as 
either  due  to  the  pressure  exerted  by  the  blade  surface  on  the  cir¬ 
cumferentially  thin  gas  stream  or  as  a  circumferentially  averaged  blade 
force  on  the  gas  stream  between  two  blades.  This  simplification  is 
first  introduced  by  Lorenz  (reference  8)  in  order  to  follow  the  flow 
along  a  given  s\irface.  The  physical  concept  involved  is  clarified  by 
Stodola  (reference  9).  Ruden  (reference  10)  further  proves,  for  incom¬ 
pressible  flow,  that  the  solution  so  obtained  will  give  an  average 
value  in  the  circumferential  direction  for  a  finite  number  of  blades, 
provided  that  the  departure  from  the  average  value  is  small.  For  turbo¬ 
machines  having  relatively  thin  blades  of  moderately  high  solidity, 
this  solution  can  indeed  be  taken  as  that  for  a  relative  stream  sur¬ 
face,  which  is  about  midway  (based  on  mass  flow)  between  two  blades. 
Because  the  circimiferential  variation  of  pressure  and  density  is  con¬ 
sidered  in  equation  (26),  it  is  better  to  refer. to  the  present  theory 
as  a  "through-flow"  theory  (following  Ruden) ,  or  "large-number-of-thin- 
blades"  theory,  instead  of  "axially  symmetric"  or  "infinite-number-of- 
blades"  theory.  (Differentiating  velocity  components  obtained  in  this 
theory  and  combining  them  according  to  equation  (16)  does  not  give  the 
true  vorticity.) 


With  this  second  condition,  equations  (24)  to  (28)  become 


1  ^(pW^r)  ^(pWg) 

r  §r  §z  ^ 


1  a  ^*f-l  -As^ 


^  ^  wj  =  0 


Si  ^u  ^  Ss  „ 

^  +  —  — N--  +  T  ^  +  F^ 


sr  -  sry  =  -  s?  +  T 

SH  „  Ss 


(29a) 


(30) 
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^-W 
Dt  “  ^ 


S(V^r) 


BW,\  ai  ,  y)  ,  T 


5H  \  I  T  t-  F 


Q  =  T  §1  =  ifwr  If  +  ¥2  If 


^  =  I?  ■"  M 


(33a) 


The  circumferential  enthalpy  or  pressure  gradient  in  equation  (3l) 
is  replaced  hy  the  term  F^.  The  corresponding  effects  in  the  other  two 
directions  are  represented  hy  F^.  and  Fg  in  equations  (  )  (  )^ 

respectively.  The  vector  F  is  perpendicular  to  W  or  related  to  the 
shape  of  the  blade  or  flow  surface  by  the  following  equations; 

FrWr  +  FuWu  +  FzWz  =  0  (34) 


Ft  dr  +  F^r  dcp  +  Fz  dz  =  0 


(34a) 


Only  six  independent  equations  exist  in  the  preceding  ^  , 

continuity  relation,  three  equations_of  motion,  one  energy  equation,  and 
one  orthogonality  relation  between  F  and  W. _  For  example,  equat 10 
(33a)  can  be  derived  from  the  equations  of  motion,  equation  of  energy 
(33),  and  equation  (34).  (See  reference  4.) 


DIEECT  AND  INVERSE  PROBLEMS 

In  a  direct  problem,  the  blade  surface  or  flow  surface  is  considered 
"bo  "be  given  *by  ‘the  equation 

S(r,0,z)  =  C  (2^) 
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and  two  relations  among  the  force  components  are  obtained  by  relating 
them  to  the  partial  derivatives  of  the  function  S: 


Ss 
S? 
1^ 
r  ^ 


(36a) 


Fu 


Ss 

1  dS 
r  'Se 


(36b) 


Equations  (29)  to  (34)  and  (36)  therefore  provide  eight  independent  rela 
tions  that  completely  define  the  direct  problem  involving  the  eight  pri¬ 
mary  variables  Wg,  Fj,,  F^^,  s,  and  I  or  H.  (The 

heat  transfer  Q  between  the  gas  and  the  blade  is  only  important  in 
the  case  of  cooled  turbine  blades  and  is  considered  to  be  given  by  cool¬ 
ing  considerations . ) 


In  the  inverse  problem,  the  blade  surface  is  to  be  determined, 
which  means  that  the  two  relations  among  the  force  components  as  given 
by  equations  (36)  are  not  available.  This  unavailability  does  not  mean, 
however,  that  the  designer  has  freedom  to  prescribe  the  variation  of 
two  variables  or  two  conditions  among  the,  variables  throughout  the  blade 
region  because,  in  order  that  the  differential  equation  (34a)  will  lead 
to  an  integral  blade  or  flow  surface  of  the  form  of  equation  (.35) ,  the 
following  necessary  (and  sufficient)  condition  of  integrabillty  must 
be  satisfied  (reference  ll) : 


f-vxf  =  0 


(37) 


which,  for  the  present  case,  reduces  to 


S  \  S  \ 
3r  ~  ^ 


(37a) 


This  condition  of  integrabillty  was  first  pointed  out  by  Bauersfeld  as 
early  as  1905  (reference  12),  but  is  neglected  in  many  recent  investi¬ 
gations.  In  effect,  it  gives  a  restriction  to  the  velocity  variation 
that  a  designer  can  specify  throu^  the  force  terms  in  the  motion 
equations  (30)  to  (32) . 


Whereas  in  the  direct  problems  two  conditions  are  obtainable  from 
the  given  surface,  in  the  inverse  problem  one  condition  on  the  surface 


16 


NAG A  TN  2302 


must  te  satisfied.  Hence,  the  designer  can  specify,  in  the  inverse 
problem,  only  one  relation  throu^out  the  hlade  region,  which  can  he  a 
variation  of  one  thermodynamic  quantity  of  gas  (Vu^  Vz;  V>  or  h),  one 
relation  on  the  hlade  surface  (untwisted  or  radial  hlade  element),  or 
one  relation  among  the  gas  properties  (constant  Mach  number),  and  so 
forth.  The  variation  prescribed  should,  of  course,  be  reasonable  so 
that  the  solution  exists. 

After  this  one  relation  is  prescribed  by  the  designer,  the  solu¬ 
tion  of  the  inverse  problem  is  quite  similar  to  the  direct  problem. 
Among  all  the  equations  to  be  satisfied,  equations  (29)  to  (34)  are 
comaon  to  both  problems.  In  addition  to  these  six  equations,  equa- 
tions  (36)  are  available  in  the  direct  problem  and  equation  (37a)  and 
the  one  relation  prescribed  by  the  designer  to  be  satisfied  are 
able  in  inverse  problem.  If,  in  the  design,  a  condition  on  the  blade 
surface  is  specified  such  as  the  blade  design  in  which  all  blade  ele¬ 
ments  are  radial,  Fj.  is  prescribed  as  zero  by  the  designer,  and  the 
integrabllity  condition  equation  (37a)  leads  to  a  simple  relation  be- 
tveen  F  ajid  F^-  These  same  two  relations  are  also  directly  given 
by  equations  (36a)  and  (36b),  respectively,  in  the  direct  problem,  and 
consequently  the  solution  for  the  two  problems  is  exactly  the  s^e. 

In  other  cases,  however,  the  solution  of  the  two  problems  is  a  little 

different . 


PRINCIPAL  EQUATION  FOR  TWO  MAIN  GROUPS  OF  DESIGN 

In  this  and  the  following  sections,  a  general  method  of  solution 
for  both  the  direct  and  the  inverse  problems  will  be  described.  From 
the  preceding  equations  the  throu^-flow  considered  herein  is  essen¬ 
tially  described  by  the  equation  of  continuity  (29)  and  the  equation 
of  motion  either  in  the  radial  direction  (equation  (30)),  or  in  the 
axial  direction  (equation  (32)).  Except  in  the  case  of  low-speed 
centrifugal  Impellers,  it  is  always  advantageous  to  use  equation  (30 ; , 
because  Fr  is  either  zero  in  high-speed  centrifugal  or  mixed-flow 
impellers  or  relatively  small  in  axial-flow  bladinp.  ^so,  _eit^r 
Sn/Sr  or  Bl/Sr  is  usually  equal  to  zero^  and  S(V^r)/dr  is  often 
given.  This  choice  is  used  hereinafter.  If  it  is  desirable  to  use 
equation  (32),  equations  can  be  developed  in  a  similar  manner. 

The  combination  of  equations  (29)  and  (30)  is  carried  out  by  the 
use  of  a  stream  function,  which  is  defined  as  follows  and  satisfies  the 
continuity  equation  (29): 
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=  rp  W2  =  rh' 


(38a) 


=  -rp*W,  =  -rb*r-h-^\ 


From  equation  (37), 


(38b) 


X 

,  *r-i  -As*  _  5^\|f 


1  1  Sh  5s* '\  5^|^ 

?  ■  72  S?  ^  5^^  5? 


•,»T-1  -As* 

^ 


1  5h  5s*^ ')  5i|; 
o2  5z  5z  /  3z 


The  succeeding  development  is  a  little  different  for  the  two  main 
groups  of  designs  to  be  considered.  In  the  first  group,  the  variation 
of  the  angular  momentum  of  the  gas  about  the  axis  of  rotation  is  pre¬ 
scribed  by  the  designer;  that  is, 

Vur  =  fi(r,z)  (4C 

is  given.  Among  this  group  are  the  free- vortex  type  in  which  fp  is 
just  a  function  of  z,  the  more  general  solid-body- rotation  type,  the 
symmetrical- velocity-diagram  type,  and  others.  In  the  second  group  of 
designs,  the  following  relation  between  tangential  and  axial  velocity 
is  prescribed  by  the  designer: 


p  =  tan  p  = 


f2(r,z) 


Among  this  group  are  the  common  blade  design  for  hi^-speed  centrifugal- 
and  mixed-flow  impellers  in  which  all  blade  elements  are  radial  with 
p  =  r  the  less  general  design  with  p  =  tan  p  =  f^{z),  which 

gives  a  practically  tintwisted  blade  and  is  most  suitable  for  cooled  tur¬ 
bine  rotor,  and  others.  The  principal  equation  will  now  be  obtained  in 
a  form  most  convenient  for  these  two  main  groups  of  designs. 
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Group  of  Designs  in  which  Equation  (40)  is  Specified 
From  equations  (9)  and  (38), 

, ,  I .  .  i  [(if)'  *  (n)' 

Differentiating  with  respect  to  r  and  z  yields 


i  =  ^[i  -  «U  (-  i  + 


■V"  (m,,  afi  -  wr  jSt 


Sh  =  ^  -  Wu  ^  -  Wm 

^  a2-W  2  5^  ^ 


u  Tj  2  Ss^ 


,-l^.T-l^*  L  -  Wr  g 


Suhstituting  equation  (43)  into  equation  (39)  and  adding  yield 


(a2-W^2)  +  (-  ^  -  i  +  «u  ^  +  a2  If-)  If  + 

oz 


|I  +  Wa  ^  +  a2  |si 
hz  Sz  dz  /  dz 


Substituting  equation  (44)  Into  equation  (30)  and  dividing  by  a^  yield 
the  following  principal  equation: 
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2  /  ar2 


+  +  |i  +  N,  $t=0 


S2  ''I  s?  ’■  ”1  sj 


where 


11  3l  „  2  ^  2  ila*  ^ 


z 


1  (  Si  ^  TJ  ,  ^2  Ss* 


With  the  variation  of 


or  W  prescribed  hy  the  designer,  the 


meridional  velocity  components  are  determined  hy  the  principal  equa¬ 
tion  (45) .  The  other  equations  are  used  to  determine  various  terms 
involved  in  the  coefficients  Lp  and  Fp.  From  the  coefficients  of 
the  second  derivatives,  the  principal  equation  (45)  is  hyperbolic  when 

the  meridional  velocity  Wjjj  =  -yW^,^  +  is  greater  than  the  speed  of 


the  meridional  velocity  Wjjj  =  y W^,  +  is  greater  than  the  speed  of 

sound,  and  elliptic  when  the  meridional  velocity  is  less  than  the  speed 
of  sound.  For  the  hyperbolic  case,  the  method  of  characteristics  can 
be  used  and  will  be  discussed  later.  For  the  elliptic  case,  it  is  con¬ 
venient  to  put  the  principal  equation  in  a  slightly  different  form,  as 
follows:  From  equation  (38), 

*  ^  S^\|f  _  1  St  _  5(loge  p*)  ^  SWj.  _  ^2^  S(loge  p*)  ^ijf 


Sr  ^^2  ”  r  5r 


Substituting  equation  (46)  into  equation  (30)  results  in 

S^t  1  M  ,  S^t  St  ^(l°Se  P*)  St  ^(l°Se  P*) 

.  2  ■  r  S?  ^  >  2  ■  [5?  S?  ^ 


rp*)^  K  S(Vur) 


+  F^  =0 
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In  this  form^  all  terms  except  the  first  three  are  taken  as  constant 
during  the  successive  improvement  of  ilr  values  throughout  the  whole 
region  in  the  numerical  solution  so  that  the  coefficients  in  the  differ¬ 
ence  equations  of  t's  will  not  vary  during  successive  cycles. 

Equation  (45)  or  (47)  is  then  the  principal  equation  for  this  group 
of  designs  to  be  used  for  a  meridional  velocity  greater  or  less  than 
the  velocity  of  sound,  respectively.  The  process  involved  in  solving 
this  principal  equation,  together  with  other  equations  in  the  inverse 
and  direct  problems ,  is  as  follows : 

In  the  inverse  problem,  or  is  given  by  equation  (40) . 

Equation  (3l)  is  first  used  to  compute  F^.  The  energy  equation  (33) 
is  then  used  to  determine  the  variation  of  s  along  the  streamline . 

The  variation  of  I  along  the  streamline  is  obtained  from  equa¬ 
tion  (33a).  Equation  (32)  is  used  to  compute  Fg  and  Fj.  is  then 
obtained  by  integrating  equation  (36a)  along  a  constant  r  line: 


(48) 


where  F  =0  at  z  =  Zq*  solution  is  then  carried  downstream  by 

equation^(45)  in  the  hyperbolic  case,  whereas  successive  sets  of 
improved  values  of  are  obtained  throughout  the  region  in  the  ellip¬ 
tic  case.  The  quantities  Wj.  and  are  then  computed  from  equa- 

^  tion  (38) . 

In  the  direct  problem,  equations  (36a)  and  (36b)  are  given.  It  is 
most  convenient  to  obtain  from  equation  (34)  as  follows: 


¥u  = 


(49) 


The  quantity  F^  is  then  computed  from  equation  (3l)  and  F2  and  F^. 
are  obtained  from  equations  (36) .  Equation  (33)  is  used  to  determine 
the  variation  of  s  along  the  streamline  and  equation  (32)  to  deter¬ 
mine  the  variation  of  I  or  H.  Equation  (30)  is  then  used  to  solve 
for  ^  as  "before. 
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Group  of  Designs  in  which  Equation  (4l)  is  Specified 

For  this  group  of  designs,  it  is  necessary  to  comhine  into 

W2  according  to  equation  (4l)  as  follows:  Substituting  equation  (41) 
into  equation  (30)  yields 


It  -  Ir "  “(7 "  i " "  I  ^  “  0  (50) 

Instead  of  equation  (42), 


t  =  I  +  -  I  r-V“ 


(51) 


Differentiating  with  respect  to  r  and  z,  combining  with  equa¬ 
tion  (39),  and  substituting  into  equation  (50)  give  the  following  fonn 
of  the  principal  equation: 


T  -u  M  -  n 


(52) 


where 


(iV) 


1  .  Ss*  1  /Si 


r  Sr  0^2  \Sr 


+  OD  r  -  W„ 


Sr; 


f\x  ^  Sq\  ^  a^-W^  /  Si 


>2  S7 


_  , ,  Ss 

2ij  2  V  Sr  Sr 


a-^W, 


+  T  ^  +  Fj.  +  ScjoWy^ 


i\T  1  4.  Sq 

^2  =  3r";[2S7^;r^"ss 


Equation  (52)  hecomes  hyperbolic  when  the  relative  velocity  is  super¬ 
sonic^  and  elliptic  when  the  relative  velocity  is  subsonic.  For  sub¬ 
sonic  velocity^  a  more  convenient  form  of  this  equation  for  computation 
is  obtained  by  substituting  equation  (46)  into  (49): 
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(iV)^ 

Sr 


1 

Sz2  '  P 


[< 


,  2\S\lf  Sp 

>sf  s? 


St  Sp 

Sz  Sz_ 


+ 


2iiCDrp*  + 


"Siir 

S7 


Sr 


Ss 
^  Sr 


(53) 


The  process  involved  in  solving  the  principal  equation  simultane¬ 
ously  with  other  equations  in  the  direct  and  inverse  problems  is  essen¬ 
tially  the  same  as  that  given  for  the  previous  group  of ^designs.  The^ 
only  difference  is  that  is  now  computed  from  equation  (4l)  and,  in 

the  elliptic  case,  is  to  be  reevaluated  after  successive  improvements 
of  t-  In.  the  design  where  radial  blade  elements  are  employed,  the  com¬ 
putation  is  considerably  shortened  with  equal  to  zero,  which  shows 

the  advantage  of  using  the  cylindrical  coordinate  system  for  this  prob¬ 
lem.  (in  the  customary  method  employing  a  coordinate  system  along  the 
streamline  and  normal  to  it,  the  blade  force  along  the  normal  is  not 
zero.)  In  this  design,  the  process  for  the  direct  and  the  inverse 
problem  is  exactly  the  same. 

The  use  of  p  =  f3(z)  in  the  design  will  lead  to  a  blade  close  to 
the  untwisted  type  if  the  hub-tip-radius  ratio  is  not  too  small.  In 
the  direct  problem  with  such  a  blade  given,  this  relation  may  be  used 
or  the  flow  may  be  more  accurately  obtained  by  using  the  equations 
given  for  the  first  group. 

The  different  character  of  the  principal  equation  for  the  two 
groups  of  designs  considered  is  interesting.  The  character  of  the  prin¬ 
cipal  equation  depends  on  the  variation  prescribed  in  the  design  or 
considered  as  given  in  the  direct  problem.  This  fact  may  be  utilized 
to  solve  some  flow  problems  in  which  the  flow  in  some  region  is  sli^t- 
ly  supersonic.  When  a  tangential- velocity  variation  is  given,  the 
equation  for  the  whole  region  may  be  elliptic. 


GENERAL  METHOD  OF  SOLVING  PRINCIPAL  EQUATION 
Elliptic  Case 

For  convenience  of  discussion,  the  fundamental  equation  (47)  or 
(53)  of  the  two  groups  of  designs  can  be  written  in  the  following  gen¬ 
eral  fcTu: 
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K  St  ^ 
r  SF  ^ 


St  Sp  St  Sp\  ,  * 


-  r^p*^M 


The  equation  is  nonlinear  even  for  incompressihle  flow.  The  equa¬ 
tion  may  he  more  conveniently  rewritten  in  a  linear  form  as 

J  ^  .  K  1+  +  ^  „  B  (5E 

::^^2  r  or 


where 


N  =  ij 


St  Sp 
Sr  Sr 


St  Sp 
Sz  Sz 


2  *2 

Sf 


and  is  evaluated  from  an  approximate  solution  at  the  start  of  the  cal¬ 
culation  and  from  the  ^  and  p  values  obtained  in  the  previous  cycle 
during  the  calculation.  For  simple  boundary  shapes  and  simple  func¬ 
tions  of  J  and  K;  it  may  be  possible  to  find  a  Green’s  fmction 
G(y^z,T];^)  with  its  proper  properties  so  that  the  solution  of  the  prob¬ 
lem  can  be  written  in  the  following  form: 


ilr(r,z)  =  dr] 


For  example,  for  flow  with  the  design  in  which  the  tangential  velocity 
is  prescribed,  the  principal  equation  takes  the  form 

5!|  .  i  ^  „  (58 

If  the  boundary  walls  are  cylindrical  surfaces,  the  total -enthalpy  and 
entropy  are  uniform,  and  the  tangential  and  the  radial  velocity  both 
vanish  at  the  inlet  and  the  exit,  the  radial  variation  of  ^  at  the 
inlet  and  the  exit  are  the  same  and  can  be  subtracted  from  "ijr,  which 
results  in 
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which  can  he  written  as 


corresponding 

in  reference  3. 
solution  for  \|f 

t  =  ti  +  r  ^  (in  (61) 


The  quantity  ( — ^7— J  is  zero  on  the  boundary  and  the 

Green’s  function  is  available  from  a  similar  eqmtion 
When  the  Green's  function  of  reference  3  is  used,  the 
is 


If  G  is  tabulated  at  several  values  of  r,  ilc-il/’q  can  be  conveniently 
obtained  by  a  numerical  double  Integration  process  on  a  punch-card 
machine.  For  curved  boundary  walls  in  the  meridional  plane,  it  is 
necessary,  in  this  method,  to  use  the  technique  of  conformal  trans¬ 
formation  to  render  the  given  boundary  shape  into  a  rectangular  one. 
Inasmuch  as  this  process  involves  a  numerical  solution  of  the  Laplace 
equation  with  the  given  boundary  shape,  it  is  found  better  to  solve 
directly  the  given  equation  (55)  with  the  given  shape  by  the  mamerlcal 
method.  Furthermore,  this  solution  will  be  the  only  choice  in  the 
general  case  where  J  and  K  are  not  equal  to  1  or  the  boundary  con¬ 
dition  is  more  general,  which  makes  the  task  of  obtaining  the  proper 
Green's  function  a  very  difficult  one,  if  not  impossible. 


In  order  to  solve  directly  the  given  equation  (55) ,  a  general 
numerical  differentiation  formula  for  first  and  second  derivatives  with 
function  value  given  at  unequally  spaced  grid  points  using  second-  and 
hi^er-order  polynomial  representation  is  required  to  give  conveniently 
and  accurately  the  finite-difference  expressions  at  the  grid  point  near 
the  cirrved  boundary,  which  is  done  as  follows:  When  the  value  of  any 
quantity  q  is  known  corresponding  to  a  n-umber  of  unequally  spaced 
values  of  the  Independent  variable  x,  denoted  by  Xq,  x,  .  .  .  x^^, 
the  variation  of  q  with  respect  to  x  is  most  conveniently  expressed 
by  a  Lagranglan  polynomial  of  the  n'*^^  degree: 


q(x) 


^  ^n+1^^)  qj  ^  0. 

^  x-x^  Ii;,+l(xj)  n 


(62) 
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where 


=  (x-Xq)(x-X2^)  .  .  .  (x-Xj^)  (63) 

^n+lC^j)  =  •  •  •  (xj-Xj-l)(xj-Xj+i)  •  •  .  (Xj-Xj,)  (64) 


2b  . 


n+1  I 


where  V  lies  between  xq  and  x^.  The  successive  derivatives  of  q 
with  respect  to  x  at  any  point  x  can  he  expressed  as  (reference  13) 

rjla  =  r  <1^,  .  ^  1  ~|  .  1r 

^  2. . /  TT I  fx - )  x-x-  2 _ /  X  XI  ^  n^ 

k/j 


1)2  ^  2*  ^  (  1  1 
U>  x-xj,  X-XJ 


+  (67) 


and  so  forth,  with 


n+1) '. 


2p  „  ^+l(^)  (n+l)/„-y  ,  ^  n:^+l(x)  _(n+2)/„^  dV 

"  “TSTTT" ‘I  (^)  +  2  -/— fv.  (u)  ^ 


n+1 


The  simmation  operation  is  very  easily  performed  when  x  is  a  grid 
point,  because  most  of  the  products  vanish.  At  these  points,  it  is 
convenient  to  write 


(D®q)x-Xi  =  +  X  (70: 

^  j=0  “ 

The  differentiation  coefficients  B  and  the  coefficients  of  the  deriv¬ 
ative  in  the  first  or  second  remainder  term  have  been  explicitly 
expressed  in  reference  13  in  terms  of  the  spacings  between  the  succes¬ 
sive  grid  points  for  general  nonuniform  spacing  throu^out  and  for  the 
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special  case  near  a  tapered' or  curved  boundary  where  oriy  the  f 
or  last  spacing  is  different  from  the  others,  using  polynomials  of 
second,  third,  and  fourth  order.  For  the  special  case,  these  coef¬ 
ficients  have  also  been  computed  for  difference 

1  29  in  intervals  of  0.01  of  the  distance  between  the  boundary  and 
the  nearest  point  and  the  other  spacings  and  are  given  “ 
ence  13.  For  spacing  lying  between  tabulated 

polation  coefficients  given  in  reference  14  can  be  used  to  obtain 
the  required  values  of  B. 

In  the  present  fluid-flow  problems,  it  is  necessary  to  cover  a 
large  region  in  order  to  reach  the  boundary  condition  that  is  always 
given  at  ^stations  far  upstream  and  downstream  of  the  ' 

order  to  reduce  the  labor  of  computation,  it  is  desirable  to  determ 
if  the  number  of  grid  points  required  for  a  given  accuracy  can  be 
reduced  ™  using  an  order  of  polynomial  higher  than  the 
seS  oSer.  A  study  of  the  remainder  terms  (reference  13)  and  actual 
experience  in  the  present  problems  show  that,  in  most  cases,  the  use 

of  fourth-order  polynomial  will  reduce  the  ^^er 

points  to  less  than  one-quarter  of  that  required  by  the  second-order 
polynomial.  In  setting  up  the  grid  pattern,  it  is  always  desirable 
to  map  the  flow  region  in  such  a  manner  that  the  distance  ® 
boundary  and  a  point  next  to  it  is  not  too  small  compar^  with  the 
other  distances,  because  the  differentiation  coefficient  becomes  ve  y 
sensitive  to  th^  small  ratio.  If  the  small  ratio  c^ot  be  avoided, 
it  is  best  not  to  include  these  points  in  the  calculation. 

When  the  grid  pattern  and  the  order  of  polynomial  representation 
have  leLSed,  the  coefficients  B  at  each  point  can  be  obtained 

from  the  table  given  in  reference  13.  Then  the  is'^Ieplaced 

tion  (55)  at  any  grid  point  whose  t  value  is  (fig-  D  is  replacea 

by  the  following  algebraic  equation: 


j=0  1 


=  0 


where  ^  and  ^  denote  the  values  of  t  along  constant  z  and 
cons'taii'b  r  lines  ^  respectively • 

The  values  of  t  along  the  hub  and  the  casing  walls  can  be  arbi¬ 
trarily  chosen,  with  the  difference  proportional  to  the  mass  flow 
-K  +  f^n  thpm  At  the  first  station  on  the  left  1-1  and  at  the  last 
staSm  2-2  in  the  right  (fig.  2),  however,  the  i  values  are 
The  boundary  condition  at  stations  1-1  and  2-2  is,  usually,  that  the 


(QOZ 


NACA  TN  2302 


27 


flow  is  parallel  to  the  bounding  walls.  When  the  hub  and  casing  walls 
are  horizontal  at  the  inlet  and  exit  of  the  machine  and  the  stations 
chosen  are  far  enough  out^  the  value  required  at  a  point  6^ 

distance  away  from  the  point  is  also  equal  to  Whether  or  not 

the  inlet  and  the  exit  stations  are  chosen  far  enough  out  will  be 
indicated  by  the  variation  of  ^lr*s  along  these  stations  obtained  in 
the  solution. 

After  equation  (Vl)  is  obtained  at  every  interior  grid  pointy  a 
number  of  methods  can  be  used  to  solve  the  set  of  simiiltaneous  alge¬ 
braic  equations.  For  hand  computation^  the  relaxation  method  developed 
by  Southwell  has  proved  to  be  superior  to  others  for  this  type  of 
equation  (references  12  and  15  to  18).  If  the  fourth-order  polynomial 
representation  is  decided^  the  calculation  can  be  carried  out  in  two 
steps  by  using  only  the  five  main  coefficients  in  early  stages  and 
then  the  residuals  are  recomputed  and  relaxed  by  all  nine  coefficients 
(reference  13) . 

If  a  high-speed^  large-scale^  digital  computing  machine  is  avail¬ 
able^  the  set  of  equations  can  be  solved  either  by  an  indirect  or 
direct  method.  In  the  indirect  method^  the  straight  iterative  method  of 
Liebmann  is  used^  wherein  is  solved  at  .each  point  from  its  sur¬ 

rounding  values  according  to  equation  (7l)  and  the  process  is  repeated 
until  the  change  at  any  point  is  no  longer  significant.  This  method 
is  simplest  to  set  up;  but  is  slowest.  A  better  method  is  to  set  the 
machine  to  perform  a  simple  relaxation  process  by  computing  the  resid¬ 
ual  at  each  point  and  relaxing  according  to  a  fixed  relation  with 
respect  to  the  amoxmt  of  residual  just  found. 

A  much  quicker  machine  method^  especially  when  a  nmber  of  solu¬ 
tions  with^  for  example^  different  inlet  Mach  numbers  are  to  be  found 
with  a  given  geometrical  shape  of  the  problem^  is  the  direct  method 
that  solves  the  set  of  simultaneous  equations  by  a  matrix  process. 

The  details  of  a  matrix  method^  which  fully  utilizes  the  great  number 
of  zerO' elements  of  the  original  coefficient  matrix^  is  given  in  refer¬ 
ence  13. 


Hyperbolic  Case 

In  the  hyperbolic  case^  the  problem  consists  in  solving  ^  from 
the  following  principal  equation^  written  in  a  general  form  for  the 
two  groups  of  designs: 
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i!i 

Bz2 


(72) 


with  the  initial  condition  that  ■>!/  and  its  normal  derivative  are  given 
on  a  curve  that  is  not  a  characteristic  line. 

From  eq,uation  (72)^  the  equation  of  the  characteristic  line  is 

2 


idz 


2K||£|  +  L  =  0 

\dz, 


The  slopes  of  the  characteristic  lines  in  the  r,z-plane  are 

^1 


/dr\ 

K 

1  [ 

Vdzyi 

J  “ 

JN 

K^-JL 


/dr\ 

-  K  +  i  [ 

\dz4 

"  J  JAI 

K^-JL 


(73) 

(74a) 

(74h) 


When  the  point  on  the  r,z-plane  moves  along  the  characteristic 
curve  I,  corresponding  to  a  small  change  dz  in  z,  the  change  in 
r  is  dr  =  dz.  Because  of  these  two  small  changes,  the  change  of 

any  quantity  q  is 


dq  =  H  dz  =  If  dz  +  ^  Xl  dz 


(75) 


or 


dq  Sq  ,  .  ^ 


(76) 


Hence,  along  2^, 


_d_  At 

dz 


s;  s)  ^  ''1  S?  155 


A.  ^  f 

dz  \Sr/  ^  \Sr/ 


,  ,  B  At 
+  ^1  5? 


+  X  S^t 
'  Sz2  1  Sz  5? 

-  4.  \ 

3z  dr  1 


(77) 

(78) 
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From  equatio 


)  and  (78), 

^  JL 

Sr  dz 

=  A  ^ 


5^i{r 

(79) 

+  A  ^  ^ 

1  5r2 

(80) 

^2.2  dz  \dz/  ”  1  dz  \5r, 


Substituting  equations  (79)  and  (80)  into  equation  (72)  yields 


d  /B' 


dz  \dz 


+  (2K-JA^)  ^  (ll)  +  (JA3_^-2K>.^+L)^ 


By  virtue  of  equations  (74a)  and  (74b),  equation  (8l)  becomes 


d 


dz 


0  ^2  ^ 


d  (M\  M  Silf  N  St 


sjy  +  j  5^  ^  j  s? 


Similarly,  along  the  second  characteristic  line 

d  /St  \  ,  d  /St\  ,  M  St  ,  N  St  _  n 
di  V5i  j  d?  teJ  j  51  j  S?  -  ° 


(82) 


(85) 


starting  from  two  points  a  and  b  a  short  distance  apart  on  the  curve 
where  the  initial  condition  is  given,  equations  (74a)  and  (74b)  give 
the  tangent  to  the  characteristic  curves  at  these  two  points  and  equa¬ 
tions  (82)  and  (83)  give  the  new  value  of  and  at  the 

point  of  intersection  c  of  the  two  tangent  lines.  The  auxiliary 
equations  corresponding  to  the  particular  problem  are  then  used  to 
determine  other  pertinent  quantities  at  the  point  c.  This  process  is 
to  be  carried  step  by  step  downstream.  The  method  is  the  same  as  for 
ordinary  two-dimensional  rotational  flow.  (For  details  of  calculation, 
see  reference  19 . ) 


APPLICATION  TO  TYPICAL  DESIGNS 

The  following  sections  include  a  brief  discussion  of  the  manner 
in  which  the  fundamental  and  auxiliary  equations  reduce  to  particular 
forms  for  several  typical  designs.  In  actual  computations,  it  is 
always  desirable  to  render  all  quantities  dimensionless.  A  convenient 
system  is  to  divide  r  or  z,  W  or  V,  p,  T,  s,  H,  or  I, 


30 


NACA  TN  2302 


'22  2 

F,  and  O)  by  r-^,  0-^,  Pt^I^  ^  ^t  At^ 

and  U-j^At^  respectively.  These  dimensionless  values  are  used  in 
the  following  equations. 


Free-Vortex  Design 


In  the  free-vortex  design,  the  variation  of  V^r  is  prescribed  as 
a  function  of  z  only.  With  a  free-vortex  flow,  the  total  enthalpy  at 
any  point  z  is  simply  related  to  the  inlet  value  by 


■q-'JK  -rr^  ^ 

Hj-Hi  - 


i  Ut" 


at  +  0?  [(vf.r* 


u 


u 


)i 


(84) 


where  (V*r*)j^  is  a  constant.  If  Hj_  is  uniform  with  respect  to  r 
and  BQ/Sr  is  zero,  SE/^r  will  be  zero  everywhere;  but  hE/dz  is  not 
equal  to  zero  in  the  rotor,  whereas  dl/Sz  is  zero  and  it  is  jt^ere- 
fore  convenient  to  use  the  system  of  equations  involving  I.  The 
principal  ec^uaiion  is  ihen 


^j,*2  r*  Sr*  Sz*2 


1  ^  ^  ^ 

Ur*-  Si^  Sz*Sz*y 


^2  ^2 
Sr* 


+  T 


Ss* 
Sr* 


=  0 


(85) 


The  auxiliary  equations  are 

*  *  _  _^  Ml 

^u^  r*p*  Sr*  Sr* 

*  wj  S(vjr*)  _^*^*  l'  St^ 

Fz  =  -  ^  Sz*'  r*2p*2  Sz* 

1  St*  ^  ^  irML  Ml  +  Ml  M1\ 

^  "  r*  Sr*  Sz*2  P*  VSr*  Sr*  Sz*  Sz*^] 


(86) 


(87) 


to 

o 

CV3 


F 


r 


r* 


a£ 


(88) 
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^  + 

(^*r*)2 

y*2 

u 

2Hf 

2Hi 

2 


[hr*)  \hz* 


r-1 


e 


-2As* 


(89) 


The  computation  can  be  started  vith  an  assumed  value  of  ijr  con¬ 
forming  to  the  boundary  shapes.  In  the  early  stage^  it  is  advantageous 
to  omit  equations  (87)  and  (88)  and  to  use  the  following  approximate 
formula  of  which  is  obtained  by  solving  the  two  equations  by 

assuming  constant  and  neglecting  small  terns: 


2  ■a(vy) 


(90) 


Density  can  be  obtained  from  the  \|r  derivatives  as  follows: 

2 

equation  (89)  is  written  as  a  =  [l  -  cp/a]  ,  where 


First, 


a  = 


*2 

I*  a)^2r*2  K_ 

^  2Hf  2Hi 


Ht 


2 

r-1 


(2r*^Ht)"^ 


CD*^r^  2 

2Hf 


y+l 

r-1 


Second,  either  op  is  computed  for  a  number  of  values  of  o  and  a  curve, 
is  plotted,  or  o  is  obtained  for  a  number  of  values  of  cp  through 
iteration  of  the  preceding  equation  and  a  table  obtained  for  equal 
intervals  in  cp.  In  either  case,  after  the  \|f  derivatives  are  obtain¬ 
ed,  cp  is  computed  and  a  or  p  is  then  obtained  from  either  the 
curve  or  the  table. 


Design  Based  on  Symmetrical  Velocity  Diagram  at  All  Radii 

Generalized  for  compressible  flow,  the  symmetrical  velocity  dia¬ 
gram  at  all  radii  is  defined  as  follows  (reference  4):  (See  fig.  2  for 
station  notation.) 
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dr^  ^  d4 


(91) 


If  it  is  desiratle  to  maintain  a  constant  total  state  along  the  radius, 
the  variation  of  V^r  with  z  is  such  that  the  same  amount  of  work  is 

done  along  all  streamlines: 


*  ^ 


(V^r  =  f(zj) 


(92) 


Hence ^ 

Sr^ 

Comhining  equations  (9l)  and  (93)  yields 

(94) 


S(vjr*)t 


=  0 


(93) 


Hence, 


(V  )e,t 


-CV^r*) 


h,t 


(95) 


(<  *)e  = 


+2 

Jl. 

2 


(^u^  )e,t~^^u^  ^h,t 


(96) 


(97) 


In  this  design,  the  quantity  Sl/5r  is  not  zero  in  the  rotor, 
hut  with  constant  work  input,  hE/br  is  zero.  Using  the  system  of 
equations  involving  H  is  therefore  convenient: 

_2^  _  1_  [M*  ^P*  +  + 

^^*2  ^2*'^  p*'\Sr*  Sr*  Sz*  Sz*  / 


*2  ^2 
r  p 


Sr* 


“ 

Vu 


S(vy) 

Sr* 


* 

Fr  + 


*  Ss* 


=  0 


*  ^ 
u 


1 


r^p* 


St*  ^(4  *) ,  St"" 

■  Sz*  Sr*  Sr*  Sz*  _ 


(98) 


(99) 
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% 


t  S(V*r*) 


Be 


Br* 


-  OD 


Bt  *  * 
*  Bz*  ^(V  ) 


Bt* 

Br* 


Br* 


*2  *2 

r  p  oz 


B^t* 

Br^2 


i_  Ml 

r*  Br* 


Ml!  _  M  /M!  ^  M!  ^ 


Bz*2 


p*  \Br^  Br* 


Bz^  Bz'* 


(100) 


*  * 
F_  =  F  r 
r  -^u 


B 

Br* 


-,* 


Ti*  * 

F^r 


d>L 


(101) 


p*2  ^ 


K*f 


iBr*; 


(M!\ 

VBz*i 


„  ^2  *2^* 
2r  p 


2 


•  2As 


(102) 


With  only  the  additional  complication  that  the  value  of  in  each 

cycle  is  determined  from  a  knowledge  of  the  streamline  in  the  previous 
cycle,  the  solution  of  this  problem  is  obtained  in  the  same  manner  as 
in  the  previous  case.  For  a  multistage  compressor,  it  is  important  to 
account  for  the  effect  of  loss  on  density  rise  by  including  the 
factor  in  equation  (l02) .  The  increase  of  entropy  can  be  estimated  by 
a  knowledge  of  the  polytropic  efficiency  (reference  4) . 


A  nonvortex-type  velocity  diagram  quite  similar  to  the  pre¬ 
ceding  one  is  the  solid-body  rotation  design  that  has  a  tangential 
velocity  varying  linearly  with  radius  in  front  of  the  rotor;  that  is, 


(v;^r^)^  =  Cr 


^2 


H 


=  2Cr> 


(103) 

(104) 


If  total  enthalpy  is  to  be  constant  along  the  radius  through  the  rotor, 
equations  (92)  and  (93)  also  apply.  Hence, 

(v!r*)j  =  +  [(v!r*)j^t  "  (v!r*)^^t] 


(105) 
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(V„r  )e  = 


(V„r 


(106) 


Comparison  of  these  equations  with  those  in. the  previous  ^esi^ 
interesting.  In  the  previous  design  the  change  of  whirl  through  the 
rotor  is  distributed  evenly  in  the  whirl  in  front  of  and  behind  e 
rotor  whereas  in  the  present  design  it  is  completely  put  into  the 
^hiS^SS  the  rotor  Except  for  this  difference,  the  calculation  is 
quite  similar  to  that  given  in  the  previous  design. 


Designs  Involving  Untwisted  Rotor  Blades 


Untwisted  rotor  blades  are  desirable  because  of  simplicity  in  mnu- 
facturing,  and  seem  to  be  the  most  practical  design  for  cooled- turbine 
THey  can  be  efficiently  need  if  the  stator  blndes  are 
deslKned  to  fit  them.  If  the  blade  Is  not  too  long,  it  can  be  desi^ 
r^Sba^L  that  Is  a  function  of  s  only.  The  principal 

equation  then  takes  the  following  form: 


(iV) 


ay 

as*2 


(l+M-  ) 


SiSilr*  Sp*  ,  ^  ^ 


Sz*  Sz-" 


2iJ£0*r*  p*  + 


r^2p*2 


(107) 


Radial-  and  Mixed-Flow  Impeller  with  All  Radial  Blade  Elements 

The  speed  of  rotation  of  the  rotor  can  be  increased  by  having  all 

blade  elements  radial.  With  F^  =  0,  the  integrabxlxty  equation  (37a) 

gives  the  result  that  F^/F^r  is  a  function  of  z  only.  Thus,  wh 

!l=r*f3(z*)  (108) 


equation  (41)  becomes 


p  =  tan  P  = 


=  -  r*  f,(s*) 


(109) 
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The  principal 


ion  is  then 


l+(r*  f  ^  ^  . 


Jhr*^ 


p*L 


Sr* 

(110) 


With  =  0^  use  of  the  three  auxiliary  equations  to  compute  is 

unnecessary. 

Equation  (llO)  is  further  simplified  with  Sl/Sr  equal  to  zero  if 
the  inlet  flow  is  of  the  free-vortex  type  or  has  no  whirl. 


Simple  Approximate  Correction  for  Blade  Thickness 

If  the  hlade  is  not  quite  thin,  it  is  desirable  to  add  a  simple 
approximate  correction  factor  b  in  the  definition  of  stream  functions 
of  equation  (38)  as  follows: 

b  ^  =  rp^W^  =  rh*^-^e"^^*W^  (ilia) 

1 

b  ^  =  -  rp*Wr  =  -  (111b) 

A  good  conception  of  this  thickness  correction  factor  can  be 
obtained  by  analyzing  the  effect  of  blade  thickness  on  the  specific  mass 
flow  along  the  mean  streamline  (based  on  mass  flow)  in  two-dimensional 
cascades.  Yet  unpublished  calculations  made  for  a  typical  subsonic  tur¬ 
bine  cascade  and  two  supersonic  compressor  cascades  show  that  the 
specific  mass  flow  pW^  on  the  mean  streamline  is  about  4  and  10  per¬ 
cent  higher  than  that  given  by  one -dimensional  calculations  correspond¬ 
ing  to  the  same  reduction  in  channel  area  for  the  subsonic  and  super¬ 
sonic  cascades,  respectively.  The  influence  of  the  blade  thickness 
also  extends  a  short  distance  upstream  and  downstream  of  the  blades. 

The  shape  of  the  mean  streamline  is  also  seen  to  follow  approximately 
the  mean  channel  line,  (but  with  lower  curvatirre) .  When  this  correc¬ 
tion  factor  b  is  used  in  equations  (ill) ,  all  the  equations  previously 
obtained  should  be  modified  by  replacing  p  by  p/b. 
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CONCLUDING  REMARKS 

Equations  of  motion  and  energy  for  unsteady,  three-dMensional  flow 
of  a  nonviscous  fluid  are  e^ressed  in  terms  of  gas  quantities  most  con¬ 
venient  for  analyzing  flow  in  turtomachines .  Entropy  change  due  t 
heat  SansL^in'^a  cLled  turbine  and  due  to  strong  shock  wave  in  super- 
sonic  flow  can  "be  taken  into  calculation. 

The  general  equations  are  simplified  according  to  the  standard^ 
assumption  for  steady  through-flow  calculation  in  turhomachines  having 
thin  blades  of  high  solidity.  The  problem  is  completely  defined  in 
the  direct  problem  with  blade  shape  given;  whereas  in  the  inverse  or 
lasiS  proMem,  »lth  tha  Inclusion  of  tha  Intagrablllty  condition  for 
ttrSairsurfLa,  either  one  flow  varlahle  or  one  relation  sccong  saver- 
al  variables  can  be  prescribed  by  the  designer. 

Throu^  the  use  of  the  stream  function,  the  continuity  equation 
and  the  equation  of  motion  in  the  radial  direction  are  combined  to 
a  principal  equation  for  the  present  problem.  The  principal  equation 
contains  some  terns  that  are  either  prescribed  or  to  be  determined  by 
other  equations  defining  the  problem.  Two  forms  of  the 
tion  are  obtained  for  the  two  main  groups  of  current  compressor 
SiS  designs  in  which  either  the  variation  of  tangential  velocity 
o^Se  varJaSon  of  the  ratio  of  tangential  to  axial  velocity  through¬ 
out  the  blade  region  are  given.  Nhen  the  tangential  velocity  is  given, 
?fe  prLciSl  equation  is^elliptic  or  hyperbolic  dep^ding  on  whether 
the  Lridional  velocity  is  subsonic  or  supersonic.  ^ 

between  the  tangential  and  axial  velocity  is 

equation  becomes  hyperbolic  when  the  relative  velocity  is  supersonic. 

A  general  method  of  solution  for  both  the  elliptic  and  the  hyper¬ 
bolic  cases  is  outlined.  Specific  applications  of  the  theory  to 
several  common  types  of  compressor  and  turbine  employing  free-vortex, 

symmetrical-velocity-diagram,  ion  fLtor 

radial -iDlade -element  designs  are  discussed.  A  simp  e 
for  blade -thickness  effect  is  also  suggested. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  October  25,  1950. 
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Figure  1,  -  Grid  pattern 


6 


rH 

CVl 


m 

U 

cu 

03 

5 

ClH 

•H 

Q 


fn 

03 


-P 


I 


o 

0) 

& 

03 

CM 

u 

P 

W 

P 

O 

rC| 

CO 

W 

1 

no 

EH 

CM 

H 

h 

O 

w 

lA 

ON 


O 


o 


Q) 

xi 

•H 

CO 

(D 

03 

Jh 

(D 

> 

03 

p:1 

a 

O 

-P 

O 

03 

-P 

CQ 


(Abstract  on  Reverse  Side) 


0 

• 

\ 

»s 

03 

0 

. 

1 

0 

1  P 

03 

0 

G 

0 

1  P 

03 

0 

P 

G 

0 

p 

1 

O  PM  P 

1  P  1 

3. 

0 

P 

03 

rH 

1  O  CM  P 

1  P  » 

p  0 

•H 

0 

CQ 

p 

P 

G 

P  P 

!T  p 

P 

P 

03  H 

p  p  a 

G 

P  p  P 

Cd  p 

P 

Cd 

5 

1 

•rl 

p  p 

0 

P  0 

0 

P 

O 

G 

p 

p  p 

1 

p  p  p 

0 

P  0 

0  P 

O 

G 

P 

o 

S 

0 

pj  03 

P 

P  tjD 

0 

o 

O  % 

0 

3  03 

P 

a  oD 

03 

o 

G 

1 

P 

P 

P 

P 

H  0  P 

0 

P 

P 

P 

G  1 

p 

Td  p 

P 

P  0  P 

0  P 

P 

P 

03 

tiD 

Pm 

P 

^  Cd 

0 

p  p  3 

P 

0 

0 

P 

03 

03  tlD 

Pi 

p  p 

0 

P  P  P 

P  0 

0 

P 

0 

•H 

P 

P 

P 

P  Sc  P 

P 

P 

O 

0 

p  p 

p  p 

P  Sc  p 

p  P 

8^ 

o 

0 

•H 

03 

*H  p 

03 

p 

P 

P 

P 

>  p 

03 

■H  p 

03 

p 

p 

P 

p 

9* 

P 

03 

•H 

P 

P 

03 

p  p 

0 

03 

P  03 

P 

P  P 

03 

p  p 

•N  G 

03 

O 

cd 

G 

0  P 

0 

p  •  0 

P 

•& 

P 

p 

O  p 

G  0  P 

0 

p  •  0 

A  A 

p 

p 

P 

G 

03 

0 

0  > 

P 

p  0 

S) 

> 

P 

o 

P  G 

03 

0  03  > 

a 

p  0 

S)  > 

P 

o 

0 

P 

0  o  > 

p 

o 

p 

H 

0 

P  P 

Si 

0  o  > 

p 

o 

p 

H 

P 

P 

•H 

•H 

>> 

G 

p  p  p 

03 

P 

p 

P  P 

P 

•H  P^ 

G 

p  p  p 

03  p 

cd 

O 

g 

P 

P 

P  P 

P 

P  P  0 

O 

G 

p 

p 

O  P 

P 

P  P  P 

•H 

P  P  0 

0  O 

G 

p 

P 

3 

*H 

O  P 

P 

^  P  P 

P 

•H 

o 

0 

P 

•H 

O  P 

P 

p  p 

P 

P 

o 

0 

G 

P 

p  p 

P 

P  P  ^ 

tc 

P 

p 

> 

G 

p  p  p 

P 

a  p 

OD 

P 

p 

> 

p 

1 

O 

P 

P  p 

P 

p 

o 

p 

0 

P  1 

O 

p  -H 

I 

P  cd  p 

P  P 

O 

p 

0 

o 

P 

P 

O 

P^  P 

0  O 

P 

p 

w 

p 

03 

o  p 

H 

O  P 

0 

0  o 

•H  »H 

CQ 

p 

0 

0 

2 

0 

P 

P  O 

P 

P  P 

p 

p 

p 

0  p 

0 

P  H  O 

P 

a  p 

P 

P 

H 

p 

p 

> 

P  03 

P 

0  H  P 

P 

p 

p 

0 

O 

p  p 

> 

P  0 

P 

O  H  P 

P  P 

3 

O 

o 

p 

o 

P 

P 

p  cd  cd 

0 

0 

03 

0 

p 

O  p 

P 

•H  cd  P 

0  0 

03 

0 

p 

G 

p 

3  'E 

P 

P  Pi  H 

P 

Pi 

>3 

G 

P  P  P 

P 

P  CMH 

P  CM 

> 

p 

P 

S  ^ 

O  P  0 

P 

0 

03 

Pi 

P 

>  P 

P 

a  ^ 

O  P  0 

P  0 

03 

P 

c 

p 

o 

cd 

P 

0 

P  P 

P 

Sop 

P 

P 

P 

o 

o 

p 

O  P 

P 

0  P  P 

P 

a  u  p 

p  p 

P 

o 

o 

G 

p 

P 

O 

p 

P  P 

O 

P 

O 

•H 

G 

p  p 

o 

H  P  P 

O 

d 

G 

•H 

Cd 

p 

p 

03 

p 

Ph 

P 

P  p  cd 

03 

>>  P 

P 

P 

p  p 

03 

P  CM 

P 

P  p  p 

03  •N 

f>3  P 

p 

p 

1 

•H 

P 

p 

P  P 

0 

G 

P 

o 

P 

P 

1  tH 

P 

P  P 

P  p 

0  G 

P 

o 

cd 

p 

P 

0 

G 

P  o 

P 

P  Pi  P 

P 

P 

P 

p 

G 

P 

A  p 

0 

G  P  O 

P 

P  CM  P 

P  -H 

«H 

p 

G 

03 

w 

P 

& 

*H 

0 

O 

^  O 

Pi  P 

G 

p 

•H 

03 

00  ti 

Pi 

P  0 

o 

O 

CM  H 

G 

p 

p 

P 

cd 

P 

>  W 

P 

O 

O 

0 

P 

• 

P 

p  p 

P 

P  >  03 

p 

l>3  P 

O 

O 

0 

1 — 1 

< 

o 

0 

P 

P  0 

G 

p  03 

03 

P 

P 

a 

Pi  P 

o 

03 

CM  -H  0 

G 

^  03 

03  p 

1 — I 

a 

CM 

p 

Pi 

P^ 

P  P 

P 

p  o  3 

P 

P 

0 

Pi 

0 

p  p 

S  p  p 

P 

•HOP 

•H  P 

0 

p 

p 

O 

P 

P 

0  P 

P 

3  ^ 

0 

> 

1 — 1 

P 

03 

P  o 

P 

P  0  P 

P 

p  P  OD 

0 

> 

p 

cd 

p 

rn 

P 

P  P 

rH 

p  p 

_ k _ 1  r,  . 

p 

03 

p 

P 

P  P 

P  P 

cd 

oJ 

0 

-p 

w 

H 

2 

0 

P 

0 

W 


o 

•H 

a 

O 

0 

4^ 


•H  (0 

•H 

U 


P  -H 


TZJ 


0  rd 
0  0 
PM  -P 
■  PI 
0 
03  03 


§ 

03 

§ 

I — i 

p 
o 

Pi  CQ 

Ph  0 
PJ 


O  -H 
G 
O 
P 


P  > 


•H 

O 

O 

P 

0 

> 


P 

PM 


0 
P! 
P  0 

H 

0  < 
P 


P 
•H  • 

G  ' 
0 


G 

03 


03 

0 

P  PI 
d  P 
0  P 
P  G 
w  cd 


p 

p 

cd 


cd 

•H 

P 
P 
0 
tiD  P 
P  cd 
cd  P 
P  P 


P  P 


cd 

P 

o 

*H 

P 

P 

P 

0 

a 


0 
p 

•H 

r— I 

PM  O 

w 

0  03 
U 


o 

p 


cd  w 

•H 

G 

P 

P  Pi 
O  O 
03  0 

P  P 

P  P 
0 

0 


p 


p  p 
o 


•H 


G 
>>  O 


P  P 


>  G 
_  ^  P 

O  P  P  P 
G  0  cd  PM 


( — i  cd 
cd 

H  p 


P 
0 
P 
0 
Cd  t3D 
p 

0  < 
p 

p 
o 


03 

03 

d 

G 

03 

Id 


G 
0 
PM  W 
m  p 


' — I 
cd 
P 


p  p  p 

0  0  P 

>  p  p 

P  P  G 


P  p 


03  P 
0  0 

P  0 

03  03 


CM 

0 

P  3 

p 

P 

P  0 

o 

p 

P 

P 

X 

cr 

o 

P 

cd  P 

•H 

G 

P 

a 

0 

P 

o 

0 

»H 

0 

P 

0 

P 

o 

0 

P 

p 

P 

OOP  p 

P 

CM 

o 

p 

P 

G 

CM 

0 

G 

P 

cd  o 

P 

0 

p 

0 

P 

CM 

A 

0 

3 

p  P 

0 

0 

0 

3 

> 

a 

cd 

P 

p 

P  P 

a 

P 

P 

p 

0 

f 

1 

•v 

03 

1 

P 

03 

0 

0 

G 

0 

H 

1 

O 

P  P 

1 

P 

1 

d. 

0 

•H 

0 

0 

1 — 1 

P 

P 

a 

G 

•H 

P 

d 

cd  p 

d 

0 

P 

ctJ 

1 

•H 

% 

p 

0 

P 

0 

0 

P 

o 

G 

P 

o 

0 

3 

03 

P 

d 

OD 

0 

O 

G 

1 

p 

P 

p 

P 

P 

0 

d 

0 

P 

P 

P 

0 

OD 

Pi 

P 

cd 

0 

Cd 

P 

cd 

P 

0 

0 

P 

0 

•H 

P 

cd 

p 

p 

P 

P 

P 

O 

0 

> 

P 

0 

•H 

p 

0 

P 

P 

d 

P 

9^ 

P 

0 

•H 

P 

p 

0 

cd 

d 

0 

0 

o 

cd 

G 

0 

0 

P 

• 

0 

d 

P 

P 

p 

p 

G 

0 

0 

0 

d 

d 

0 

S) 

> 

p 

O 

0 

P 

P 

rS 

0 

o 

> 

•H 

o 

P 

P 

p 

P 

•H 

•H 

a 

P 

p 

P 

0 

d 

0 

o 

P 

P 

P 

P 

P 

p 

P 

p 

0 

0 

o 

G 

d 

P 

P 

•H 

o 

*rH 

P 

cd 

P 

P 

•H 

o 

0 

>3 

G 

P 

p 

p 

4-^ 

d 

d 

OD 

d 

•rH 

> 

p 

1 

O 

P 

•rH 

1 

•H 

p 

d 

d 

d 

o 

p 

0 

o 

P 

P 

C 

^  f — ! 

0 

0 

o 

-H 

•rH 

03 

0 

0 

P 

P 

o 

P 

d 

p 

P 

f — \ 

P 

P 

> 

“H 

0 

cd 

O 

p 

p 

P 

d 

d 

Cj 

o 

P 

C’ 

P 

r-M 

*H 

cd 

cd 

0 

0 

0 

cb 

p 

>3 

G 

p 

cd 

P 

P 

P  P 

P 

P 

> 

P 

P 

R 

O 

0 

-1 

0 

00 

P 

d 

o 

P 

P 

0 

•H 

.r-t 

P 

a 

G 

p 

P 

P 

•H 

o 

o 

p 

1 — i 

P 

O 

H 

P 

P 

o 

d 

O 

*H 

G 

p 

P 

0 

P 

CM 

p 

p 

p 

0 

03 

•s 

>3 

P 

P 

Cd 

1 

•H 

P 

P 

P 

d 

p 

0 

G 

43 

o 

0 

P 

P 

0 

O 

P 

o 

p 

cd 

p 

P 

P 

•H 

•H 

P 

G 

P 

P 

"M 

0 

o 

o 

P  P 

O 

p 

•H 

0 

P 

0 

p 

1 — i 

> 

03 

4-^ 

>3 

c6 

o 

o 

0 

P 

• 

P 

o 

0 

p 

•H 

0 

G 

p 

03 

0 

p 

p 

a 

P  P 

< 

p 

P 

CM 

P 

P 

cd 

•H 

O 

0 

•H 

p 

0 

P 

0 

P 

O 

P 

p 

0 

Cd 

p 

d 

P 

OD 

0 

> 

( — 1 

0 

0 

P 

P 

P 

P 

d 

d 

>3 

9^ 

0 

0 

0 

cd 

C3 

P 

d 

•H 

*H 

P 

P 

■0 

P 

G 

d 

>3 

0 

•H 

03 

o 

4^ 

O 

•  rH 

p 

> 

0 

.H 

G 

P 

P 

G 

•H 

d 

p 

d 

P 

fj 

'H 

d 

P 

0 

P 

•H 

•H 

>3 

•H 

p 

o 

0 

>3 

o 

p 

P 

0 

P 

•rH 

0 

P 

P 

P 

P 

P 

G 

V 

d 

P 

rH 

o 

0 

OD 

G 

P 

P 

G 

o 

o 

P 

P 

0 

p 

•rl 

0 

r-H 

0 

0 

cd 

0 

0 

0 

P 

P 

p 

G 

> 

G 

0 

< 

P 

0 

p 

P 

p 

OD 

P 

o 

a 

O 

•H 

P 

cn 

•H 

1 — 1 

O 

P 

p 

a 

O 

o 

P 

rH 

P 

P 

P 

0 

a 

•H 

G 

0 

G 

0 

0 

P 

M 

P 

P 

4) 

0 

> 

d 

> 

•H 

•rH 

o 

l>3 

0 

0 

3 

P 

P 

1 — i 

cd 

0 

o 

0 

J — 1 

fp 

p 

d 

P 

p 

P 

EH 

P 

P 

P 

♦rH 

P 

1 — ! 

0 

p 

p 

0 

0 

•H 

0 

O 

p 

Cd 

cd 

0 

> 

d 

P 

00 

p 

> 

P 

00 

a 

P 

d 

•H 

•H 

G 

p 

• 

CM 

0 

p 

H 

d 

P 

d 

0 

o 

P 

P 

0 

<! 

03 

p 

d 

>< 

P 

o 

d 

0 

P 

•H 

G 

4^ 

a 

p 

0 

0 

0 

•H 

o 

0 

•H 

0 

P 

0 

d 

o 

p 

CM  P 

0 

P 

p 

P 

OD  1 — 1 

d 

•ri 

P 

o 

p 

p 

P 

P 

P 

G 

p 

0 

G 

d 

0 

o 

P 

0 

p 

p 

P 

0 

0 

p 

P 

0 

0 

•rH 

•H 

0 

0 

03 

3 

p 

0 

OQ 

> 

S 

cd 

EH 

P 

4J 

P 

P 

a 

P 

•H 

p 

